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In this paper we consider a novel statistical inverse problem on 
the Poincare, or Lobachevsky, upper (complex) half plane. Here the 
Riemannian structure is hyperbolic and a transitive group action 
comes from the space of 2 x 2 real matrices of determinant one via 
Mobius transformations. Our approach is based on a deconvolution 
technique which relies on the Helgason-Fourier calculus adapted to 
this hyperbolic space. This gives a minimax nonparametric density es- 
timator of a hyperbolic density that is corrupted by a random Mobius 
transform. A motivation for this work comes from the reconstruction 
of impedances of capacitors where the above scenario on the Poincare 
plane exactly describes the physical system that is of statistical in- 
terest. 

1. Introduction. The recovery of objects, for example, densities and func- 
tional thereof, based on noisy indirect observations, otherwise known as 
statistical inverse problems (see, e.g., [23]), is scientifically of intense inter- 
est. The literature is vast and we mention only a few selected papers. Most 
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of the work is concerned with deconvolution on Euchdean spaces. Promi- 
nent approaches are based on wavelet and wavelet-vagulette expansions (see 
[1, 24] and [39]) or, on singular value decompositions (see [5, 26] and [33]), 
where for the latter block thresholding techniques lead to adaptive estima- 
tors (see [8]). Minimax rates in deconvolution have been investigated by Fan 
[15], and many others. Recently, oracle inequalities have been proved in [9]. 
The specific problem of boxcar deconvolution and its link to Diophantine ap- 
proximation have been investigated in [21] and [22]. Other methods include 
the linear functional strategy (see [16] and variants thereof). Very popular 
are Fourier series estimators which have been known for a long time, and 
often they become particularly simple because on Euclidean spaces they can 
be treated with kernel methods (see [15, 19] and [44]). We note that even 
though our approach also utilizes Fourier methods on groups, due to the 
hyperbolic geometry, the resulting estimator cannot be treated by kernel 
methods thus complicating our endeavor considerably. 

In this paper we provide a novel methodology for statistical object recov- 
ery on the Poincare upper half plane which we call the problem of Mobius 
deconvolution. Here the group of Mobius transformations is given by all 
fractions of the form (az + b)/ {cz + d) for a complex number z and a 2 x 2 
matrix with real entries a,b,c,d of determinant 1. The metric which is in- 
variant under these transformations is the hyperbolic metric (to be specified 
in the next section) which will replace the Euclidean metric in a natural way. 
In fact, a key observation is that this problem can be tackled by generalized 
Fourier methods similar in spirit to the Euclidean case. The development of 
the theory builds on the foundations laid out in [43], Chapter 3. However, 
extending Euclidean arguments to this manifold is challenging since the hy- 
perbolic space is noncompact, the (hyperbolic) geometry is non-Euclidean 
and the group of Mobius transformations is noncommutative. A fundamen- 
tal technical difficulty comes from a lack of a dilation property that does 
not extend over from the Euclidean case. Despite these difficulties, a generic 
element is the Riemannian structure on a manifold, on which a Laplacian 
can then be defined. This together with the generic Euclidean approach as 
outlined in [33] will be the foundation to what will be presented below. 

In addition to the theoretical interest of this novel scenario of Mobius 
deconvolution, there are important practical applications as well. One par- 
ticular situation occurs in alternating current circuit analysis and design 
whenever signals travel through circuit elements as well as whenever geome- 
tries of waveguides change. A simple example of the latter is a connector 
to a coaxial cable, say. Here, the so-called "reflections" are modeled on the 
complex unit disk and the corresponding "impedances" occur in a complex 
half plane. Usually, these reflections or impedances are not directly visible 
but observed through other electrical devices, such as a "two-port" which in 
turn is modeled by Mobius transformations. In particular, a class of so-called 
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"lossless" two-ports can be identified with 2x2 matrices of determinant 1. 
One particular aspect of Mobius deconvolution is related to the temporal de- 
cay of impedances of capacitors whereby the above scenario on the Poincare 
plane exactly describes the physical system that is of statistical interest. 
Other applications include the field of electrical impedance spectroscopy, as 
well as electrical impedance tomography. In the former, measuring varying 
impedances due to variable ion transport through biological membranes is 
currently of high interest in view of pharmaceutical drug design (see [14] as 
well as [41]). In the latter, in a noninvasive and radiation-free way, medical 
imaging can be cost effectively accomplished by measuring skin-impedances 
(cf. [6]). Indeed, for successful reconstruction, control of various errors is of 
paramount importance (see [17]). There is also work in higher-dimensional 
hyperbolic spaces with respect to medical imaging (see, e.g., [29]). 

In statistics there is also some recognition of the Poincare plane and its hy- 
perbolic geometry particularly so because the parameter space of the Gaus- 
sian distribution (with unknown mean and standard deviation) is this space. 
Furthermore, it has been shown by several authors that the Riemannian met- 
ric derived from the Fisher information is exactly hyperbolic (see [25, 32] for 
details). Obviously, location and dispersion parameters of arbitrary distri- 
butions and random estimators thereof can be viewed within the Poincare 
plane. Curiously here, the family of Cauchy distributions play a specific role 
as being equivariant under Mobius transformations (see [34, 35] for this and 
its consequences for parameter estimation). Based on the above, techniques 
from the hyperbolic geometry of the Poincare plane are developed exclu- 
sively from a parametric point of view (cf. [36]). As far as the authors are 
aware, our contribution is the first attempt at nonparametric developments. 

We now summarize the paper. Section 2 is a preliminary section which 
introduces the notation along with the Helgason-Fourier analysis needed for 
this paper. Following this. Section 3 presents the main results. In Section 
4 we focus on computational aspects of Mobius deconvolution illustrating 
the ideas through simulations. To this end we introduce in addition to the 
hyperbolic Gauss the hyperbolic Laplace distribution. In Section 5 we go 
into explicit detail with respect to the Mobius deconvolution problem for 
statistically recovering the temporal decay of impedances of capacitors as 
outlined two paragraphs above. We will briefly sketch the background; how- 
ever, if the reader is well versed in this field, then one can start from Section 
5.4 where we examine a data set that was acquired through collaboration 
with the University of Applied Sciences (Fulda, Germany) that depicts the 
physical system of this paper. In particular, we are able to identify random 
impedances when only their impedances viewed through random capacitive 
two-ports are given. Following this, technical details of the Poincare upper 
half plane and the proofs of the main theorems are collected in Appendices 
A and B. 
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As usual for two function g and /, write / x (7 if f{x) = 0{g{x)) and 
g{x) = 0{f{x)) for 2; — )• 00 or x — )• 0, depending on context. 

2. Preliminaries. In the following let R and C denote the real and com- 
plex numbers, respectively. Furthermore, the group of real 2x2 matrices of 
determinant one is denoted by 

(2.1) SL(2,M) := 1^^= ■.a,b,c,deM.,ad-bc=l 

This defines the group of Mohius transformations Mg : C — t- C by setting for 
each 5GSL(2,M), 

(2.2, 

where MgMh = Mgh for 5, ft, G SL(2, M). Let 

(2.3) ]H:={zGC:lm(z) >0} 

be the upper half plane where "Re(z)" and "lm(z)" denote the real and 
imaginary parts of a complex number z, respectively. Then for each g G 
SL(2,M), the Mobius transformation Mg is a bijective selfmap of H. More- 
over, for arbitrary z, z' G EI there exists a (in general not unique) g G SL(2, M) 
such that z' = Mg{z'). 

The action of SL(2,R) on H, which is rather involved, is further discussed 
in Appendix A. It will be used in the proof of the lower bound in Appendix 
B.3. For the following we note that Mobius transformations preserve the 
family of vertical lines and circles centered at the real axis (cf. Figure 1). 
This is a consequence of the fact, that Mobius transformations leave the 
cross ratio 

{Zl - Wl){z2 - W2) 

c{Zi,Z2.,Wi,W2) = 7 TT r 

[Zi - Z2){Wi -W2) 

invariant. For a detailed introduction (cf. Nevanlinna and Paatero [38], 
Chapter 3). 

The deconvolution, or statistical inverse problem of reconstructing the 
density of a random object X on H, of which we only see a version Y cor- 
rupted by an independent random error e on §L(2, M) can now be formulated 
as 

(2.4) Y = MAX). 

A natural geometry for (2.4) is the given by the hyperbolic distance on H 



l + \l\c{z,z',z',z) 

a{z, z ) = log 



1 - ^J\c{z,z',z',z)\ 
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since for this distance the space of isometries of H is precisely the group 
of Mobius transformations, meaning that d{Mg(z),Mg{z')) = d{z,z') for all 
g G SL(2,R) and z, z' € H. Here, 1 = x — iy denotes the complex conjugate 
of z = X + iy. The corresponding hyperbolic measure is chosen such that: 

(i) its area element dz agrees with the area element dx dy of Lebesgue 
measure at z = i, and 

(ii) it is invariant under Mobius transformations. 

In consequence, the Radon-Nikodym derivative of the hyperbolic area ele- 
ment with respect to the Lebesgue area element at w = u + iv is given by 
v~'^ which is the determinant of the Jacobian 



of a Mobius transformation M yielding M(i) = (ai + b)/{ci + d) = u + iv. This 
can be verified with the complex derivative M'{i) = {ci + d)^ = Ux + ivx and 
the Cauchy differential equations Uy = —Vx,Vy = Ux] cf. Terras [43], Chapter 



In addition, in order to properly define below in (2.6) a convolution of a 
density on H with a density on §L(2,M), a compatible bi-invariant Haar 
measure dg on §L(2,R) is chosen in Appendix A. 

Hence, X and Y are random complex numbers in the upper half plane H 
equipped with the hyperbolic geometry, and e is a random isometric self-map 
of EI applied to X by (2.2). The problem of the Mobius deconvolution can 
be made precise as follows. A density on the upper complex half plane with 
respect to the hyperbolic measure is called a hyperbolic density. Densities 
on SL(2,M) are taken with respect to the Haar measure dg. 




III. 



At z = X + iy & M we have hence the hyperbolic area element 



(2.5) 



dxdy 
dz := — 5—. 

yZ 



y 




X 



Fig. 1. The shortest connection (a geodesic segment) between two points G H in the 
hyperbolic geometry is either a vertical line segment [ifKe{z) —Ke{z')J or an arc on the 
circle through z and z with center on the real axis. 
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Problem 2.1. Under the model (2.4) estimate nonparametrically, the 
hyperbohc density fx of X from the hyperbohc density fy of Y when the 
density fe on SL(2,]R) is known. 

We note that this setup assumes underlying i.i.d. Xi,...,Xn corrupted 
by i.i.d. errors ei, . . . , En, also independent of Xj (j = 1, . . . , n) giving obser- 
vations Yj = M^^ {Xj ),j = 1, . . . ,n. 

We will base our work using Fourier or singular value decomposition meth- 
ods that are common for the Euclidean case (see [5, 26] and [33]) and the 
fact that the densities of (2.4) are related by the convolution 

(2.6) fy{z) = {f,*fx){z)=: [ fe{g)fx{Mg-.{z))dg. 

JSL(2,R) 

This fact is a consequence of (A. 4) in Appendix A. 

From here on, we will make the abbreviation SL(2) := SL(2,M), as well 
as to write (2.2) as simply g[z) or gz for g S §L(2) and 2 G H whenever the 
context is clear. 

2.1. Fourier analysis on the Poincare plane. For purposes of Mobius 
deconvolution for Problem 2.1 we sketch the Helgason formulation of hy- 
perbolic Fourier calculus which can be found in more detail in [43], Chapter 
3.2. The Helgason-Fourier transform of / G C^(BI) with the latter being the 
space of real valued functions with compact support in M with derivatives 
of all orders, is defined as the function 

nf{s,k):= [ f{z){lMk{z))fdz 

analytic for (s, fc) G C x SO(2) where overline denotes complex conjugation. 
Here, 

k = ku=( ""^"^^ ^'''''^ GSO(2) CSL(2) 
y — smii cosu J 

is naturally identified with u G [0, 27r) acting on H as the Mobius transfor- 
mation 

, ^ , , z cos u + sin u 

Mk{z) = : 

cos u — z sm u 

as defined in (2.2) (cf. Appendix A). Note that for all s G C, z = x + iy ^ 
and z — )■ {lm(k(z))Y are eigenfunctions with corresponding eigenvalues 
s{s — 1) of the Laplace-Beltrami operator 
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on H. With the spectral measure 

dr = — ^ttanhfvrt) dtdu 
on M X SO(2) the inverse Helgason-Fourier transform is given by 

(2.8) f{z)= [ r~'\f(}- + it,K){lm{k{z)))'/^-^''dT, 

where = —1. The following result justifies these definitions: mapping to 
the Helgason-Fourier transform extends to an isometry L'^{M,dz) — )■ x 
SO{2), dr); that is, we have the Plancherel identity 



(2.9) / \f{z)\'dz 



u=0 



^/( ^ +it,ku\ dr, 



where we denote the space of square integrable functions over some space 
by L^. We note that / G L'^{n,dz) is SO(2)-invariant if and only if Hf G 
L^(M X §0(2),dr) is SO(2)-invariant. Thus, for numerical computations 
([43], pages 141 and 149) for an §0(2)-invariant function /, transforms and 
inverse transforms can be considerably simplified 

"W/Q+^i^ =27ry" /(e~''f)P_i/2+it(coshr)sinhr(ir, 

(2.10) f{e-'^i) = l-j ^/Q+ft^P_i/2+,t(coshr)ttanh(7rt)dt 
with the Legendre function 

Paic) ■.= ^ {c+^C^-lC0s{(t>)Tdct>. 

Throughout this work we will use the following assumptions: 
(D.l) all densities are square- integrable 

fxjY^L\M,dz), /,GL2(SL(2),d(7); 
(D.2) the error density is bi-invariant 

f,{agh) = f{g) G SL(2), a, 6 G SO(2); 

(D.3) fx G Fe,{Q) for a Sobolev bah 

-^a(Q) = {/ G L\n, dz) : ||A"/2/f < Q} 
with a > 1 and Q > 0. 
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Here, A'^/^/ denotes the unique function h € iF'iW^dz) with T-Lh{s,k) = 

As detailed in Appendix A, the isometry SL(2)/SO(2) -;>EI:5SO(2) ^ 
Mg{i) preserves the action of SL(2). Hence a density satisfying (D.2) can 
be regarded as an §0(2)-invariant mapping EI — )■ IR. In particular, in case of 
(D.l) and (D.2), the Helgason-Fourier transform T-Lf^{z) is well defined, we 
have 

(2.11) nfY{s,k)=n{fe * fx)(.s,k)=nfe{s)-nfx{s,ky, 

[43], page 149. One final assumption to be made is the following. 
(D.4) 3 constants /3,7,Ci,C2 > such that 

Ciexp|-^| < |^/e(s)| <C72exp|-^| Vs = ^ + ft, teR. 

As an example, the hyperbolic Gaussian-distribution (see Section 3 below) 
satisfies (D.4). 

Of course, any density on the upper half plane (or on the unit disk) 
can be rescaled with respect to hyperbolic measure. One example, using the 
normalized squared absolute cross ratio c{z, 6,z, 6), has been kindly provided 
by one of the referees, 

1 f\z-z\\e -6\y , Aa'^dxdy 

vr V \z-6\^ ) 7r((x - /i)^ + (y + cr)^)^ 

with a hyperbolic parameter = /i + icr G H; for example, one could take 
X and y as suitable estimators of the location and dispersion parameters 
of another distribution (cf. [34, 35]). Note that this density is not SL(2)- 
invariant; rather it is equivariant with respect to the SL(2) action on both 
variable z and parameter Q. In the context of this research one is interested 
also in SL(2)-invariant densities. Such can be generated from suitable densi- 
ties on y G [l,oo). Moreover, SL(2)-invariant functions additionally fulfilling 
(D.3) can be obtained by applying the inverse Helgason transform (2.10) to 
suitable functions on s = ^ + zt, — oo<t<oo. If the function is even in t, 
then the inverse Helgason transform thus obtained is real. It is, however, 
not necessarily nonnegative. As a consequence of P_i/2+ji (coshr) > for 
t = 0, nonnegativity can be obtained if the function tends sufficiently fast to 
zero as t — t- oo. Numerical experiments indicate that one may consider for 
T > —1/4 and a > 1 a suitable multiple of a power of a Cauchy density in 
the spectral domain 

1 1 

(2.12) nh^^ris,k)(x 



(t-s(s-1))° (T + l/4 + t2)<^ 
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giving an invariant hyperbolic Laplace density ha,T G J~a{Q)- This density 
can be lifted as in Appendix A giving a bi-invariant density ha^r on SL(2). 
In particular ha^^n * ^02, t2 ^ J^ai+a2{Q) for ti,T2 > —1/4:, ai,a2 > 1 and 
suitable Q >0. 

We will not elaborate further on this topic, but we mention that motivated 
by our research and by many potential applications, the task of generalizing 
non-Gaussian distributions to hyperbolic spaces may lead to a new field of 
challenging research. 

3. Main results. Let us begin with the definition of the Helgason-Fourier 
transform of the generalized derivative of the empirical distribution /^^(z) = 
^ X]j=i ^Yj i^) where Yi,. . . ,Yn is a random sample in H 



1 " 

(3.1) nf^'\s,k) = -Y,{MHYjW. 

Obviously 

(3.2) Enfir\s,k) = nfY{s,k), 

where "E" denotes expectation. We estimate the Helgason transform of an 
§0(2)-invariant density as well by an SO(2)-invariant estimator 

(3.3) nfPis) ■.= nfP{s,i) = -^im(y,)^ 

with the identity element I S §0(2). 

From the Helgason-Fourier transform (3.1) we build an estimator by using 
(2.11) and the inverse Helgason-Fourier transformation (2.8) with a suitable 
cutoff T > 

(3.4) fp^\z) := I r" ^^IU^{lHku{z))f"^^' dr 

J\t\<TJu=0 nfe{l/2 + lt) 

for the density fx- This is well defined if Tifs 7^ is bounded from below 
on compact sets which is guaranteed under assumption (D.4). Even though 
we consider in this section the general case, we note in view of (3.3) and 
(2.10) that for the estimation of an SO(2)-invariant density fx, we can use 
the simpler 

As the first main result we have: 
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Theorem 3.1. For fx,fY and fs satisfying (D.1)-(D.3), and Tifs ^ 
bounded from below on compact sets, there is a constant C > not depending 
on T, a, Q and n such that 

nfx^^-fxf<C sup 

\t\<T 

as n — )• oo . 

If the corruption by error is smooth enough, or equivalently if the asymp- 
totic rate of the decay of its Helgason-Fourier transform is suitable, the 
cutoff T can be adjusted appropriately to obtain the following rates. 

Theorem 3.2. Suppose that fx.fy and f^ satisfy (D.1)-(D.4). Then 
by letting T = logn — ^ log(logn))"'^/^ where i] > 2(a + l)//3 

m^x'^'^ - /xf <Qi^llognj (1 + o(l)) 
as oo where a is from condition (D.3). 

The optimal rate of a power of logn in case of error smoothness (D.4) 
is in agreement with Euclidean results. On the real line, condition (D.4) 
corresponds to supersmooth errors for which Fan [15] estabhshes the same 
type of rate. This rate has also been established by Butucea and Tsybakov [7] 
in case of additionally supersmooth signals. For a scenario corresponding to 
our setup on compact Lie groups, see [28] and more general on any compact 
manifold, see [27], where similar rates have been found. 

The above results are minimax in the sense that the rate of convergence 
is matched by a corresponding lower bound. We have the following theorem. 

Theorem 3.3. Suppose that fx.fy and f^ satisfy (D.1)~(D.4). Then 
for some constant C > 0, we have 

inf supEll/" - fxf > C(logn)-2°/^ 

as oo, where the infimum is taken over all estimators /" and the supre- 
mum over all fx G J~aiQ)- 

Recall the Gaussian density gp on the real line with zero mean and vari- 
ance 2p > can be characterized as yielding the solution of the heat equation 



nfe{- + it 



— + QT- 
n 



2a 



(A - dp)u = 
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with initial condition u{z,0) = f{z) by 

(3.5) u{z,p) = {gp* f){z). 

Similarly on M, the density gp giving the solution of the heat equation by 
(3.5) is also called the Gaussian density for M. Here, the Laplace-Beltrami 
operator A would be defined by (2.7), and the convolution in (3.5) would 
be defined as in (2.6). Using (2.11) for SO(2)-invariant / and u it is easily 

seen that 'Hgp{s) cx e'^^'^"^^^. Consequently, in terms of assumption (D.4), the 
Gaussian density satisfies (3 = 2 and 7 = We have the following result. 

Corollary 3.4. For fx and fy satisfying (D.1)-(D.3) consider cor- 
ruption according to a Gaussian distribution fe = 9p- Then by letting = 
^ [log n — r] log(log n)] where rj>l + a, 

E||/.?''^^-/xf >^(logn)"° 
as n —7- 00 gives the optimal rate of convergence. 

4. Computations and simulations. In this section we elaborate on com- 
putational aspects, simulations and, in particular, discuss the Gaussian dis- 
tribution on H. We begin by first discussing methods for choosing the trun- 
cation parameter. 

4.1. Estimating truncation parameter. A popular technique for data- 
driven choice of a truncation parameter is least squares cross-validation 
(see [13] or [45], Chapter 3.3). We will discuss how that technique can be 
adapted to our setting. For a given random sample Yi , . . . , y„ , an optimal 
cutoff T = T* > minimizes the mean integrated squared error 

T„* = argminjEf / {fp^\z)f dz\ - 2¥.( f fx{z)f'~^^^\z)dz 

Instead of deriving a minimizer of the above we content ourselves with min- 
imizing a suitable estimator. Obviously, ^^{f'^^''^\z)Y dz is an unbiased es- 
timator of the first term. Let 

and therefore choose 

T. := argmin { [ {fp'^kz))' dz - ^± f^^^^A , 
T>o \Jm J 

which is an estimate for an optimal T = T*. 
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Alternatively, we can use the result of Corollary 3.2 and set 

1/^ 



T- 



-logn- - log (log n)'' 



We are aware of the fact that cross-validation in general suffers from too 
large variability and, of course, more involved parameter selection methods 
could be generalized here as well (see, e.g., [11, 12, 37, 40] and [42] among 
many others). However, we do not pursue this issue any further in this paper. 

4.2. Simulation of the Gaussian distribution. For simulation we use the 
analog gp of the Gaussian distribution on the upper half plane introduced 
above. Recall that by a more subtle argument (see [43], pages 153 and 155), 
the inverse transform is obtained in polar coordinates [for any k G SO(2)] 

gp{k{e i))=gp{e z) = — =3 V2e p' / , . , = 

yJA.T^p Jr vcosho — coshr 

2tt sinh r ' 

that is, for a px-Gaussian distributed SO(2)-invariant random object X on 
in and an §0(2)-invariant subset ^ C H, 



F{XgA}= [ 9pxir)dr 
Jw>0:e-'-ieAnm\ 



{r>0: e-'-jeAnH} 

[see (A. 2)]. Hence, in order to simulate X from an invariant Gaussian dis- 
tribution we simulate rx ^9px ^ "^x uniform on [0, 27r); then 

X = kux° Rrx (^) = ^ux (e"''^i). 

Note that gp^. G J^a{Q) for all q > with suitable Q = Qa > 0. That is, fx 
satisfies condition (D.3). 

Similarly, in order to simulate e from a bi-invariant /Og-Gaussian distribu- 
tion on SL(2) we consider ~ gp^ and ku^ , ku'^ independent and uniform on 
[0,27r).Then 

Y = ku,o R^^ o ku'^{X) 

[see (2.4) and (A.3)]. 

According to Theorem 3.1, (3.1) and (3.4) we can then estimate the den- 
sity fx by 
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100 independent hyperbolic 0.1-gauss 
distributed deviates 



• A** ' * 



Each deviate on the l.h.s corrupted by an independent 
bi-invariant 0.05-gauss Moebius deviate 



1000 independent hyperbolic 0.1-gauss 
distributed deviates 



Each deviate on the i.h.s corrupted by an independent 
bi-invariant O.Q5-gauss IMoebius deviate 




10000 Independent hyperbolic 0.1-gauss 
distributed deviates 



Each deviate on the l.h.s corrupted by an independent 
bi-invariant 0.05-gauss Moebius deviate 



Fig. 2. Three simulated data samples of n — 100 (top row), n = 1000 (middle row) and 
n — 10,000 (bottom row) independent random data points on the upper half plane. Left row: 
original independent invariant px -Gaussian distributed data points, right row: transformed 
data points under independent n bi-invariant p^-Gaussian distributed §L(2) transforma- 
tions (right). 



By SO(2)-invariance it suffices to estimate for 
estimate ^^-^ (r) by 



e "^i only, hence we 



fx (^) 



27rsinhr/j^''^^(e 



with the integral simplified as in (2.10). 

In the following simulation we consider an original distribution with px = 
0.1 under a corrupting Mobius transformation distributed with = 0.05. 
From this we create three data sets with different sample sizes: n = 100, 1000 
and 10,000. Figure 2 shows the original ^i, . . . ,Xn and the corrupted data 
Yi = (-'^i), • • • , = Mfr^ {Xn) in cartesian coordinates in the upper half 
plane. In Figure 3 we show the corresponding densities times hyperbolic 
area on [0,oo). Note that these are then densities in the usual sense; that 
is, their integrals with respect to Lebesgue measure on [0,oo) are 1. The 
density estimation by deconvolution has been obtained from the observed 
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Densities for 100 independent triais Deconvolution with) cutoff T = 3.39 




0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 



Densities for 1000 independent trials Deconvolution witti cutoff T = 4.43 




Densities for 10000 independent trials Deconvolution with) cutoff T = 5.51 




Fig. 3. Left-hand side: population and empirical densities times hyperbolic area [corre- 
sponding to Qp-^ (r) and Qpx+pei''^)] along the first polar coordinate r of the data depicted 
and described in Figure 2. Top row: n — 100, middle row: n — 1000 and bottom row: 
n = 10,000. Right-hand side: additionally the respective estimate times hyperbolic area: 
p^'"^^(r) of the original density by Mobius deconvolution. The corresponding optimal cut- 
off parameters T = Tn have been estimated by least squares cross-validation as in Section 

4.1. 



data Yi,. . . ,Yn by the proposed method. For the deconvolution, since only 
the optimal rate 

. ^ \i/4 f 2.19, (n = 100), 
r« — logn = < 2.42, (n = 1000), 
^ ^ {2.61, (n = 10,000), 

is guaranteed by Corollary 3.4, we have used the estimate via least squares 
cross-validation as proposed in Section 4.1. 

4.3. Simulation of the hyperbolic Laplace distribution. Using formula (2.10) 
directly with (2.12) to simulate Laplace (a, r)-deviates for Laplace dis- 
tributed SL(2) error corruption, we obtain results similar to the ones re- 
ported above. Due to the oscillation of the Legendre polynomials to be 
evaluated, however, the computational time is much longer. In analogy to 
Theorem 3.2 we have the upper bound 0(n~"/^^^"^) for the choice T = 

^l/(2(a+l))_ 
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(a) 



(b) 



Fig. 4. Basic circuit models for signal processing. Left: one-port, right: two-port, (a) 
Serial circuit with generated voltage ua, generator impedance Zc and load impedance Zl; 
(b) circuit of Figure 4(a) with a two-port inserted between generator and load, depicting 
input (a.) and output waves (h.) at the two-port and at load. 



5. Impedance density estimation in AC driven circuits. For the conve- 
nience of the reader, we begin this section with a review of classical electri- 
cal engineering theory specifically tailored to the application of hyperbolic 
statistics in mind. For the underlying engineering terminology we refer to 
standard textbooks such as [10]. More mathematical approaches are ex- 
plained in [2, 20] and [43], Chapter 3. In the following we rephrase this 
problem in the language of statistics. We are then able to identify a typical 
problem as a novel inverse problem in hyperbolic space. 

Here, general Mobius transformations appear with complex coefficients a, 
6, c, d in (2.1). Moreover, hyperbolic space materializes in the form of the 
upper half plane H, the open unit disk D := {u; G C : |tt;| < 1}, and the open 
right half plane -iU := {Q G C:Re(C) > 0}. With the notation of (2.3), all 
are related to one another by Mobius transformations, the first is usually 
called the Cayley transform 

w=C{z):= — — , z = i- , tC = z, 

Z + l 1 — w 

w = - , C = • 

C+1 l-w 

5.1. Complex impedance in AC circuits. We begin our discussion with a 
one-port, a single load impedance serially inserted in a circuit of a voltage 
generator and its impedance [see Figure 4(a)]. Recall that voltages, currents 
and impedances in an alternating current (AC) circuit are modeled by com- 
plex numbers; otherwise, a loss of alternating real voltage u:t^ UQCOs{ujt) 
over a load giving a phase shifted current j :t ^ jQCOs{ujt -\- (f)) would re- 
sult into an awkward time dependent real resistance u{t)/j{t). In complex 
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notation, the ratio of voltage u{t) = uqc^^^ over current j{t) = j^e^^'^^^'^^ is 
constant and called impedance 

Z:=^ = ^e^^eC. 
jit) jo 

Its real part is called resistance, the imaginary part is the reactance. For 
example, under an AC-voltage uq cos(a;)t, a serial circuit of a resistor with 
direct current DC-resistance R and an ideal capacitor with capacitance C 
features an inverse impedance (called admittance) of = + icoC. In 
fact, in realistic scenarios, the resistance is positive, thus Z G —iM, and the 
boundary (the imaginary axis) corresponds to ideal (lossless) impedances. 



5.2. Reflections and characteristic impedance. We now assume that our 
circuit features a generator generating the open circuit voltage uq with in- 
ternal impedance Zq and a load with impedance Z^ as depicted in Figure 
4(a) with total impedance Z = Zq + Zl according to Kirchhoff 's circuit law. 

Inspired by the wave model, voltage loss ul over and current flow Jl along 
the load is considered to be the superimposition of an incoming (denoted by 
"+") and a reflected wave (denoted by "— ") in such a way that each single 
wave satisfies Ohm's law with a common characteristic impedance Zc. Since 
the reflected wave propagates into a direction opposite to the incoming wave, 
we have the ansatz 



UL = u\ + u^ and jL=j+-j^ 



with Ohm's law 



Uf Ur uf+Uj 

-T = Zc = -^ and i =Zl. 

3l Jl 3l 3l 

The specific decomposition or equivalently the choice of Zc is arbitrary in 
many applications, and it will be guided by imposing additional conditions. 
Usually Zc is taken positive, or at least chosen such that the normalized 
impedances Z := Z/Zc will be again of positive real part, that is, Z G —iM. 
The analog for the right half plane of the Cayley transform, (5.1), yields 
then reflection coefficient of the load 

(5.2) r, := ^ = 5 = = ^ = C{-^Z,) 



as an element of the unit-disk D. Of course, there is no reflection if Z^ = Zc- 
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5.3. The chain matrix. We are now in a position to investigate the 
generic scenario of signal transmission through a two-port [see Figure 4(b)]. 
Among others due to linearity of the Maxwell equations, voltages and cur- 
rents {jl = ai, j2 = 02) on either side of the two-port have a linear relation- 
ship governed by a so-called impedance matrix Z 



Z 



32 ) \ Z2I Z22 ) \h) \U2 



ch 



For given circuit parameters, the coefficients of the impedance matrix can 
be easily computed. For example, Zn = ^|j2=o is the well-known input- 
impedance (by inserting a load of infinite impedance the right-hand side 
becomes an open circuit with j2 = 0). In most applications it turns out that 
Z is symmetric; the corresponding two-port is then called reciprocal. 

One easily verifies that the chain matrix which is usually denoted by ch 
(the Russian letter "cha") relating (^2,-^2) with (iii,ji) is given by 

n2 \ J_ / Zii det(Z)\ / U2 \ ^ ( ui 

-32 J ^21 V 1 ^22 ) V -32 ) V 

[in contrast to the mathematical literature, the engineering literature tends 
to use a transmission matrix relating [112,32) with (ni,ji) instead, with 
reversed j2 = —02 in Figure 4(b)]. An advantage of the chain matrix over the 
impedance matrix is that the former is well defined for the limit Z21 — ?• 00, 
for example, for an ideal coil in series with the load. 

Again, only using lossless (i.e., purely imaginary) impedances (such as 
ideal inductances, transformers and capacitors) guarantees that the corre- 
sponding chain matrix has real diagonal coefficients and imaginary coeffici- 
ents elsewhere. Moreover, for cascaded two-ports (i.e., several two-ports in 
serial connection), the resulting chain matrix is just the product of the indi- 
vidual chain matrices. By linear algebraic decomposition of §L(2,]R) it can 
be shown that every lossless two-port of lumped elements can be modeled 
by cascading combinations of two-ports involving only inductances, trans- 
formers and capacitors (see [20], page 18). Note that j2 = 02 = — = —ji 
in Figure 4(b). Hence ch defines a Mobius transformation relating load im- 
pedance with the input impedance of the two-port 

fry \ chiiZL + chi2 chiiU2 - Chi2j2 ui 

IVl(.h[^L) — = = — = ^1- 

ch2lZL + ch22 Ch2in2 - Ch22j2 jl 

Here, Zi = ui/ji is the impedance of the load Zl = —U2/j2 as viewed 
through the two-port. As a consequence we make the following remark. 

Remark 5.1. Serial cascading of lossless two-ports is equivalent to the 
action of the Mobius group SL(2, M) on the i-fold iZ^ G H of load impedances 
Zl. 
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We are thus led to the statistical inverse problem (cf. Problem 2.1). 

Problem 5.2. Estimate the load impedance Zi^ when only the impe- 
dance Z{^x) viewed through the two-port can be observed where e = Xo Mch o 
X^^ € SL(2) is assumed to be of known distribution. Here, X : C — ^ C : z — )• 
denotes the multiplication with i. 

In^conclusion we note that one may as well consider normalized impedan- 
ces Z = Z/Zc or equivalently reflection coefficients, (5.2). Then the mapping 
for the normalized impedances goes as follows: 

Zi _ diuZL/Zc + chi2/Zc 

Zc ZcChi2ZL/ Zc + Ch22 ' 

5.4. Estimating resistances seen through electrolyte capacitors. It is well 
known that over the duration of years properties of electronic equipment 
change due to wear-out effects of various elements. In particular, electrolyte 
capacitors have a tendency to loose capacitance. In effect, older electronic 
devices deviate from original calibration and may feature nondesired side- 
effects; for example, field strengths of transmitters may grow stronger than 
tolerated. 

In an application of our method we consider a series of n = 150 measure- 
ments of random resistors of 15 $7 provided with an accuracy of 10 percent by 
the manufacturer (they range from 13.5 ft to 17.7 il) viewed through 30 ran- 
dom capacitors at 1 kHz taken at the Department of Electrical Engineering, 
University of Applied Sciences, Fulda, Germany. These originally identical 
lossy 22 iiF capacitors have been collected from over ten year old electro- 
nic gear. For the impedance measurements the LCR-Bridge "HM8118" has 
been used that comes with an accuracy of 0.3% guaranteed by its producer 
HAMEG. We model the ?-fold of the impedance Z of these capacitors with a 
hyperbolic Gaussian-distribution at unit impedance, that is, Zi = ZcMp^[i) 
with a suitable characteristic impedance Zc and a random hyperbolic Gaus- 
sian p-distributed Mobius transformation Mp^ . Measurement of the capaci- 
tors gives 0.0004 corresponding to a spread of roughly 4.8%. Our goal 
lies in the reconstruction of the one-dimensional resistances R solely from 
the observations 

W= ' =^ 
l/R + l/Z Z + R 

and the known dispersion of the corruption as posed in Problem 5.2. To 
this end we apply Mobius deconvolution to the model 
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Random resistors seen 
through lossy capacitors 



°8 oo 



hyperbolic mesh 




1 ] ] ] I ] 

2.5 3.0 3.5 4.0 4.5 5.0 5.5 



deconvolution along 
fixed angles of mesh 




random resistors seen 
through mean capacitor 



0.2 0.4 06 0.8 1.0 

distance to center 



Fig. 5. Random resistors seen through random capacitors. Left top: original measure- 
ments. Right top: hyperbolic polar mesh points at which the Mobius deconvolution was 
computed. The bottom left image depicts the deconvolved densities times hyperbolic mea- 
sure along fixed angles. The red and blue line goes along the red and blue mesh points. For 
verification, the one- dimensional distribution of the resistors seen through the single mean 
capacitor is depicted in the bottom right image. 



with suitable characteristic impedances Wc and Re- For this apphcation the 
Euchdean means have been chosen as characteristic impedances. Alterna- 
tively, a better approach may be to use hyperbolic intrinsic means (see [3] 
and [4]). Figure 5 shows the observations W in the left top corner. Mobius 
deconvolution is computed along the hyperbolic polar mesh-points depicted 
in the right top corner (Figure 5). Below in the bottom left corner the de- 
convolved densities along fixed angles of the mesh are depicted. The angle 
depicted in red shows highest density followed by the angle depicted in blue. 
The location of the two dominating directions depicted with the same colors 
in the right top corner (Figure 5) is in high agreement with the location of 
the impedances of the resistors seen through the mean capacitor depicted 
in the bottom right image of Figure 5. Indeed, one can say that with the 
few measurements available, we were able to reconstruct the nature of the 
unobserved elements, namely resistors with impedances distributed along a 
one-dimensional subset in the complex half plane. 



20 



HUCKEMANN, KIM, KOO AND MUNK 




Fig. 6. Simulation ofn— 1000 random resistors seen through n random capacitors. Top 
left: observed impedances, characteristic impedance denoted by "x." Top right: unobserved 
underlying normalized resistors seen through unobserved mean capacitor. Bottom left: ra- 
dial density times hyperbolic measure of nonobserved resistors. Bottom right: deconvolved 
densities times hyperbolic measure along fixed angles of mesh in Figure 5. 



In Figure 6 the above scenario is more prominently reenacted in a si- 
mulation of n = 1000 measurements using Rc = min(i?) and Wc from the 
preceding example (depicted by "rr" ) . We show observed measurements and 
unobserved capacitors in the top row as well as original density and the deco- 
nvolved densities along respective angles as in Figure 5 in the bottom row. 
Obviously the distribution of the unobserved resistors is quite reasonably 
recovered along the grid. 

APPENDIX A: POLAR COORDINATES AND CONVOLUTION 

In this appendix we focus on the right and left action (they are different 
due to noncommutativity) of the special orthogonal group §0(2) on §L(2) 
giving rise to polar coordinates and to M viewed as the quotient with respect 
to one of the actions. In fact, the action of §L(2) on H can be naturally 
viewed in polar coordinates (cf. Figure 7) of which we will make extensive 
use in the proof of Theorem 3.3 in Appendix B.3. 
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Fig. 7. The polar coordinates {r,u) of z £M are obtained from the unique circle hyperbo- 
lically centered at i (star) containing z; that is, this circle is symmetric to the imaginary 
axis and intersects it at points of form e~^i and e^i (x-crosses). Rotating z by the hyper- 
bolic angle —u along this circle, the lower point e~^i is obtained. The polar coordinates 
(u2,R,u\) of g £ §L(2) rotate z by the hyperbolic angle ui along the above circle to obtain 
zi, rescale z\ by e^^ to obtain zi = e^^z\ and subsequently rotate 22 by the hyperbolic 
angle U2 along the unique hyperbolic circle through Z2 hyperbolically centered at i yielding 



In contrast to a metric on a manifold in the usual topological sense, a 
Riemannian metric is a metric in every tangent space, that varies smoothly 
with the offset of the tangent space. Thereby, every Riemannian metric 
defines a metric on the manifold in the usual sense and a unique volume 
element, called the Riemannian volume, giving rise to a unique measure on 
the manifold. For more details we refer to Lee [31], Chapter 3. 

If the Riemannian metric on SL(2) underlying the Haar measure dg is 
chosen such that the natural Riemannian quotient metric on H yields the 
hyperbolic measure dz on H, then we are able to lift densities on the bottom 
space IHI to the top space SL(2) to obtain (A. 4) yielding (2.6). 

Let us begin with the observation that the hyperbolic measure (2.5) is 
the Riemannian volume element of the extension of the standard Euclidean 
metric in the tangent space of EI at z = i in a left-invariant way under the 
action of SL(2) on H. Similarly, we equip SL(2) with the Riemannian metric 
obtained from the left SL(2)-invariant extension of the standard Euclidean 
metric in the tangent space of the unit matrix I S SL(2). We denote the 
corresponding Riemannian volume element which defines a left-invariant 
Haar measure by dg. According to [43], Exercise 19, page 149, it is also 
right invariant, that is, 



As mentioned before, for arbitrary z,z' GM there exists a g £ §L(2) such that 
z' = Mg{z). Given one such g, any other g' £ SL(2) satisfies Mg'{z) = z' if, 
and only if, g~^g' G SO(2). In particular, Mg{i) = i if, and only if, g £ §0(2). 



5(^)- 
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This entails that the fohowing mapping of the quotient space SL(2)/§0(2) 
due to the right action of SO(2) is well defined and bijective 

SL(2)/SO(2)^BI \ 
^ ' gm{2) ^ Mg{i) j ■ 

Since the mapping preserves the action of SL(2), the natural Riemannian 
quotient metric of SL(2)/SO(2) is isometric with the the hyperbolic metric 
of M. 

Next, consider the left action of SO(2) on SL(2). This projects to a left- 
action on H giving rise to polar coordinates u G [0, 2tt), called the hyperbolic 
angle and r > (cf. Figure 7) of 

, , / cosn smu\ ( e"''/^ \ 

Z=Mh O Md(i) with fc„ = . ],iir={ r/9 • 

^"^^ ' \^-smn cosny ' \^ e'/^y 

For ease of notation, for the entire paper we identify M^.^ with ku and Mji^ 
with Rr, respectively, such that kuoRri}) = ku{e~'^i). Since for every z € EI\ 
{i} we have that k(j,{z) = z if, and only \i^ = Omodvr, polar coordinates cover 
the hyperbolic plane twice. As a consequence, z G EI has polar coordinates 
with r > uniquely determined and u unique modulo vr if r > 0. Thus, the 
hyperbolic area element (2.5) transforms to 

(A. 2) dz = sinh(r) dr du. 

Polar coordinates can also be defined on SL(2): every element g S SL(2) has 
a decomposition 

(A. 3) g = kuRrku', r G [0, oo), n, n' G [0, 27r), 

with uniquely determined r > 0; if r > then u and u' are also uniquely 
determined modulo vr (in fact, one of the two is unique modulo 27r). In 
view of the isometry (A.l), this gives our choice of Haar measure in polar 
coordinates 

dg = sinh(r) dr du du' . 

If g,g' are independent random elements in SL(2) with densities /i,/2 
continuous with respect to Haar measure, we have for the probability that 
the product is contained in a measurable subset A C SL(2) by left-invariance 
of the measure that f{gg' £ A) = /^(/i * f2){o)da with the convolution of 
/i and /2 given by 

(/i*/2)(a):= / h[g)f2{9'^a)dg. 

J§L(2) 

In general, convolutions over noncommutative groups are noncommutative. 

Suppose now that Z is a random quantity on HI with density /2 conti- 
nuous with respect to the hyperbolic measure. Using polar coordinates, this 
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density lifts to a right SO(2)-invariant density /2 on SL(2) f2{kuRrk'y) '■= 
f2{kuRr{i))- Hence, convolutions of a density fi on S'L(2) with a density 
/2 on H can be well defined by lifting to a right §0(2)-invariant density on 
§L(2) 

(A.4) (/i*/2)(z):= / Mg)f2ig-'a)dg 

with any a G SL(2) giving Ma{i) = z. This convolution is commutative if 
either fi is bi-invariant or if /2 is SO(2)-invariant. 



APPENDIX B: PROOFS 

B.l. Upper bound: Proof of Theorem 3.1. In order to measure the per- 
formance of fx''^'^ we consider the mean integrated squared error 

nfx - = nfx - E/i"'^¥ + \mp^^ - /xf 

with the usual variance-bias decomposition. The assertion of Theorem 3.1 
then follows from the following more detailed lemma. 



Lemma B.l. For fx^fv o,nd satisfying (D.l) and (D.2), and T-Lfe 
bounded from below on compact sets, there is a constant C > independent 
of T and n such that 

Hfx ^fx II ^inf|,|<^|^/,(l/2 + .t)|2 n- 
If, additionally, fx satisfies (D.3), then 

l|E/?''^^-/xf <QT-2-. 

Proof. We first note that by (2.8), definition (3.4), since the right-hand 
side is in L^(EI, dz), 



•'^ ' / nfe{l/2 + it,k) ^ 



Here I denotes the indicator function. Hence by the Fubini-Tonelli theorem, 
(3.4), (3.2), (2.8) and (2.11) 

-H{Efp^^)(^^+it,k 

^ , ™-2-E1^/«(l/2 + it',fc,, 



\t'\<T Ju'=0 

nfe{i/2 + it') 
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(B.l) =nfx(^l+it,ky^^r,T){t). 



Deduce from (3.1), 

E 



(B.2) 



1 



nfri - + it,k 



+ 



Im(EA:(y)) - |?^/y(l/2 + it, k)(' 



n 



In addition, lm{Kk{Y)) = H/y(l,A;) implies that 



(B.3) 



\lm{EkuiY))\du: 



\nfY{i,ku)\du<c 



with a suitable constant C > 0, since Tify is analytic. Thus, using the Plan- 
cherel identity (2.9), (2.11), the Fubini-Tonelh theorem, (3.2), (B.l), (B.2) 
and (B.3) [by hypothesis (D.l) and (2.9), = p/yf < oo], we have 

indeed for the variance 



E /l/ 



in,T) 



E/ 



(n,T)|2 



X 



dz 



t<\T\ J«G[0,27r) 
1 1 



\nfe 



dT 



< 



noir Jt<\T\Jue[o,2Tr) 
C 



lm{Eku{Y))-\nfYil/2 + it)\' 

\nfe{i/2 + itr 



ttanh(7rt) dtdu 



inf|,|<r|H/e(l/2 + it)|2 n 

with a constant C > involving neither n nor T. 

In the next step we similarly estimate the squared bias under the addi- 
tional assumption (D.3) using also (2.9) and (B.l) 

\Efp^^ - fx\^ dz 



^ / / \nfx?dT 

J\t\>T JMg[0,27r) 

u=2n 



1 

8^ 



t\>T Ju=0 



- + it 
2 



•tit 

2 



MOBIUS DECONVOLUTION 



25 



T-(-fx ( 2 + ^" 



ttanh(7rt) dtdu 



<QT 



'2a 



□ 



B.2. Optimal rate: Proof of Theorem 3.2. If /g satisfies (D.4) we have 
the upper bound 



(B.4) log C + log Ci + 2 log T - log n + ^T^ (^1 + 



/3/2 



for the logarithm of the variance term (cf. Lemma B.l). A sufficient condition 
for convergence of the variance term while T = T(n) — )• oo, is that (B.4) tends 
to — oo. Hence, T is of form 



T{n) = ( - logn- -A{n] 



with A{n) —7- +00 at a rate lower than that of logn. A short computation 
gives the rate 



log(logn) — A{n) 



for (B.4). In case of optimality this rate must be larger or equal to the 
logarithmic rate of the upper bound of the bias term in Lemma B.l which 
is then 

log (logn). 

In consequence the rate of A(n) is r/ log(log n) with r/>2(l + a)//3as asserted 
by Theorem 3.2. 



B.3. Lower bound properties: Proof of Theorem 3.3. Recall the decom- 
position in polar coordinates from Appendix A (cf. Figure 7). The general 
idea of proof goes as follows. Define the dilation of H by 

(B.5) H\k{e-'i)) ■.= H{k{e-^'i)f-^^^^P^i,2{cosh{5r)). 

smhr ' 

For arbitrary SO(2)-invariant H and 5 > 0, is obviously also SO(2)- 
invariant 

To derive a lower bound for estimating fx in Lp' norm, we follow a classical 
scheme which has been condensed in [15], pages 1261 and 1262 (cf. also [18], 
pages 1555 and 1556). The adaption of this scheme to the Poincare plane, 
however, is not at all obvious and will be the subject of the following sections. 
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After some elaborate preparation in the following two sections we take a pair 

/O e ^a{Q), fn G J^a{Q), for wllich 

fn = f0 + CH6-''H\ 

where S = 6n (cf. Sections B.3.3 and B.3.4 below). Then, in Section B.3.5 we 
show that 5 can be chosen such that 

poo 

(B.6) x'(/. * /O, /. * fn) ■■= / (/. *fo-fe* fnfUe * /o)"' dz < -. 

Jo n 
In consequence by (3.3) of [15], there is di > such that for any estimator 

r of /x, 

(B.7) sup P/{||r-/x||2>||/0-/n||2/2}>dl, 

/xG{/0:/n} 

which gives with (3.4) of [15] a lower bound 

sup i^lir - fxf > ^U/O - /nll2 >^ S-'^WH'h- 

The choice of 6 in (B.17) at the end of Section B.3.5 in conjunction with 
jjff'^]]^ X WHW^ from Lemma B.5 then yields the rate 

asserted by Theorem 3.3. 

B.3.1. Convolution equation in polar coordinates. 

Lemma B.2. Suppose f is bi-invariant and h is SO{2) -invariant. Write 
z = kuRrii) , g = kaRskjj ^§L{2) and dg = shih s dads d/3 . Then, 

{f*h){z) = 2TT / f{e-H)h{e'k.^{e-''i))smhsdq)ds. 

J^=0 Js=0 

Proof. Since k^^ku = k^^a observe that 
{f*h){z)= [ f{g)h{g~'z)dg 



f{kaRski3)h{k_f3R_sk-akuRr{i)) 

a,/3G[0,27r) J s=0 

X sinh s dads df3 
/ f{e~^i)h{R^sk^^Rr (i)) sinh sdcpds. 

1=0 Js=0 



□ 
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Lemma B.3. Define rj{r, s, 4>) and R{r, s, 4>) by k^(r,s,4>)^~^^'^'^'^^i = e^k^^ x 
e~^i, where < r]{r,s,(j)) < 2-k and R{r,s,4>) > 0. Suppose (j) € [0,27r) and 
r, s > 0. Then, 

|r — s| < i?(r, s, (/)) <r + s. 

Proof. Let tp = 20. From [43], page 125, we take 

^ , — sin -0 sinh r + i 

cos ij: sinh r + cosh r 

Let k^(j.^g^^^e~^^'^'^''^\ = X + iy. Then using [43], page 150, 

cosh(i?(r, s, 0)) = 

2 y 

1/ _ sin^ 0^ sinh^ r + 1 

e (cos-^smhr + cosnrj + e 



2 \ cos 'ijj sinh r + cosh r 

Set t = cos-i/). Since 

sin^ sinh^ r + 1 (1 — t^) sinh^ r + 1 
cos ■0 sinh r + cosh r t sinh r + cosh r 
= cosh r — t sinh r, 

we have 

cosh(ii(r, s, (/))) = ^{e~^(coshr + tsinhr) + e*(coshr — t sinhr)} 
= — tsinhrsinhs + coshrcoshs. 

Since —1 < t < 1, cosh(i?(r, s, 0)) has the maximum cosh(r + s) at t = —1 
and the minimum cosh(r — s) at t = 1. Suppose that cosh a < cosh 6 with 
6 > 0. This imphes that 6 > |a| since 

b> a, if a > 0, 

b > —a, if a < 0. 

The desired result follows from this and the assumption R{r,s,(p) > 0. □ 

The following lemma is an immediate consequence of Lemma B.3. 

Lemma B.4. Suppose H is SO{2) -invariant and H{e~^i) is monotoni- 
cally decreasing in r. Then, 

H{e'k_^e~'^i) < H{e~\^-'\i) fors>0,(/)e [0,27r),r > 0. 
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B.3.2. Dilation. 



Lemma B.5. Suppose 
(B.8) 'HH{^ + it) = /ort^[i,l]- 

Then for 5 — )• oo, 



(B.9) 



1 



1 



1 t 



and for f^ satisfying (D.4) with 7 = 1, with a constant C > 0, 

I \fe*H'\^<Ce-'('/'^''. 

Proof. Let's start with an alternate representation of the Legendre 
function from [43], page 158, and a specific derivative 



^-i/2+ii(coshr) 



^(coshr) 



\/2 cos(tn) du 

vr Jq yj cosh r — cosh u ' 

a2 



[dtY 

V2 



l/2+it 



(coshr) 



V? du 



TT Jq \/ cosh r — cosh 



u 



In the next step we make use of [43], Exercise 28(b), page 158. For fixed 
t, K > 0, and large (5, that is, r = k/(5 — t- 



i/2+it( cosh - 1 - -^ol 



i=o 



(i!)2 V2<^y 



Similarly, 



-P-l/2+it/5(cOshK) = P_i/2(cOshK) ^ ( - | + 
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UH^ (^ + it^ = H{e-^i)P,y2 (cosh + sinh(K) dn 

/•CO 

/ H{e~'^i)P^i/2+it/5{coshK)smh{K) dhi 
Jo 



2tt 
T 



Moreover since 



t tanh(7rt) dt 



1 

47r 



1 



nH[ - + iu 



■tanh(7r(5M) 
tanh(7r«) 



iitanh(7ru) du 



and taking additionally into account that 6 is large 



|A°/2i^5||2 



t tanh(7rt) dt 



oo / 1 

^ +4 



■tanh(7r(5u) 
tanh(7rii) 



utanh(7ru) du 



<C(52"||A"/2/i-||2. 
Since by hypothesis 

\nfe{l/2 + it)\< C2 exp[-|l/2 + itf] = C2e-(^/^+*')'^', 

we obtain 

(B.IO) sup \'HUl/2 + i6u)\<C2e-^^l^+''l^^"' <C2e-^'l^^\ 

mG[1/2,1] 

It follows from (B.8), (B.9) and (B.IO) that there is a constant C > such 
that 

j \fe*H'\'' = j \nH' {1/2 +it)\''\nfe{l/2+it)\''dT 

5~^\'HH {1/2 + iu)\^\'Hfe{l/2 + i5u)\^{5u) tanh(7r(5n)5 du 

'1/2 

This completes the proof of Lemma B.5. □ 
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B.3.3. Bound on go = fe * fo- Choose /o as 

fo{k{e-'i)) = fo{e-'i) = ^(coshr)"'^ for a > 1. 



Lemma B.6. fo is SO (2) -invariant and 

f oo /■27r 

fo{ku{e~^ i)) sinhr dr du = 1. 

r=0 Ju=0 

Proof. By the definition of fo, it is SO(2)-invariant and 

foo /•2tt 

fo{kui&~^i)) sinhr dr du 

r=0 Ju=0 

POO 

= (a — 1) / (coshr)~"sinhr (ir 

Jr=0 

POD 

= (a-l) / x~"'dx = l. 

Jx=l 

Lemma B.7. Suppose that f^ is bi-invariant and 



□ 



/ fs{e ''i)smhsds > - 
Jo ^ 



for some positive constant C^- Then, go = fe* fo satisfies 

( 2^(a - l)2'^-le-2'^^ for r > C, 
\2^(a- l)2'^-ie-2«^% forr<Q 



goie "'i) > 



Proof. Note that < s < Cg < r implies r + s <2r whereas < s, r < 
Cg, implies r + s < 2Ce- It follows from these facts, > 0, Lemmas B.2 and 
B.3 that 

/.*/o(e-^i) 

/'27r /"OO 

= 27r/ / /,(e-'^i)/o(A:,(,,,,^)e-«(^'^''^)i)sinhsds# 
Jo Jo 

= {a-l) / /£(e~^i)cosh""(i2(r,s,0))sinhs(isd(/) 
Jo Jo 

>27r(a-l) / fe{e ''i) cosh "(r + s) sinhsds 
Jo 



> 27r — ^— ^ > 2^(a - l)2'^-ie-2«- 
- 2cosh"(2T) - ^ ' 
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with r = r for r > and t = Cs for r <C£. □ 

B.3.4. Bound on * . We note from [30], page 188, 

P„i/2(coshr) = ^^ir(tanh(r/2)), 

' 7rcosn(r/2j 

where the complete ehiptic integral of the first kind is defined by 

fvr/2 

K{t) 



\/l — sin^ (p 
In consequence there is a C > such that for all r > 
(B.ll) P_i/2(coshr) < C. 

Define 

/i5(e-''i) = 5e-('"o-i)'5^ for r > 0. 

When mo > 1 and 5>0, ||;U5||oo<<5 and ns{e~^'i) is monotonically decreasing 
in r. 

Lemma B.8. Suppose \H{e~^i)\ < Ce"'"^'' with niQ> 1 and 6>l. Then, 
there is a constant C > such that H^{e~'^i) < Cfis{e^^i). 

Proof. By Taylor's expansion, 

(e-)'5_(e--)'5 = 5^<5-l(gr_g-r) 

for e-*^ < r/ < Since rj^"^ < e^-^"^)^ < e''^ 



sinhr 

Hence with C > from (B.ll), 



\H'{e-''i)\ = |F(e-'5^i)|^^^^i^P_i/2(cosh(<5r)) 

smh r ' 



smh r 

Lemma B.9. Suppose that 

/ fe{e~H) smh sds < Ce-^^-^o)"" for < < I and C > I + Co, 

J r-^or 

H{e~'''i) is bounded and monotonely decreasing in r and satisfies the tail 
condition 

(B.12) H{e-''i)<Ce-"'°'' /or mo^o > 

Then, for 5 >1 and r > 0, 

{fe*H^){e-'i)<C6e-^^-^°>. 
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Proof. Set TZi = {s:\r- s\ < ^^r} and TZ2 = {s :\r - s\ > (,or}. It follows 
from Lemmas B.2, B.4 and B.8 that 

l-if,*H')ie-^■^) 

Z27V /"OO 
/ /e(e-^i)i/^(fe^(,,,,^)e-^('^'^'<^)i)sinhs#ds 
^1=0 Js=0 

/'27r /"CO 
J^=0 Js=0 

(B.13) <2ttC J fe{e~H)fis{e-^''~'^i)smlisds. 

Set 

Ij = / fs{e~'^i)fj,sie~^^~'^k)smh.sds for j = 1,2. 
Jn, 

Since mo > ^/^q > 1, 

||/^5||oo < S. 

Observe that 

h<\\lJ's\\oo fe{e~H)smhsds 

J \r—s\<^or 

(B.14) < C7(5e~(«~«»)^ 

From the tail condition (B.12), 

l2 = [ /e(e-"i)((5e-('"o-^)^l^-"l)sinhsds 
< ^g-{mo-i)55or f f,{e-H)smhsds 

J \r-s\>£,or 

(B.15) < 5e-('"°-^)''^o^ 



Since r > 0, mo^o > C and 5 is large, 



e-(«-?o)r- 



Combining (B.13), (B.14), (B.15), we have the desired result. □ 
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B.3.5. Chi-square distance: Proof of (B.6). With the above notation, 
choose a pair of densities 

/o(A:(e-^i)) = /o(e-'-i) = ^(coshr)-'^ and fn = fo + ChS-'^H' , 

where a satisfies 1 < a < 2, and H satisfies the hypotheses of Lemmas B.5 
and B.9. By choosing Ch close to 0, we have /o,/n S J^aiQ) for ah large 6. 
Let go = fe* fo and gn = fe* fn with satisfying the hypotheses of Lemma 
B.9 and (D.4) with 7 = 1. 

The distance between go and gn is defined by 

xigo,gn)-= dz = 2TrCHS / -. — — smhrdr 

J go J go{e "^i) 

and with a suitable constant M > guaranteed by Lemma B.7, by Lemmas 
B.5, B.8 and B.9 

— '^^ sinh r dr 

go{e ''t) 

<(f +/ )'^-f-7''%inh.dr 

< — — — {fe*H (e ^i)) sinhr dr + — / — sinh rdr 

where C = 27r(a - 1)2""^ For the last equahty we set M = e^l'^^^l'^^^ , //q = 
2(^ — ^0) — 2a — 1 and [X\ = min(^o/2, (2 — a)). Then indeed 

M2a^-2iS/2f ^ ^~{2-a)i5/2f ^ o{e~^^-'') and d^Nr^^ = 0{e~^^-"'). 
Hence, 

(B.16) x'(go,gn) = 0(5-2"e-^^W2)^). 

Letting e~^^^^^'^^^ =n~^, or equivalently 
(B.17) 5 = 2/i7^/'^(logn)i/^, 

we conclude that the right-hand side of (B.16) is of order o(n~^), that is, 
(B.6) is proven. 
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